Quantum interference in dirty ci-wave superconductors 
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The local differential tunneling conductance on a Zn impurity in a disordered d-wave supercon- 
ductors is studied. Quantum interference between many impurities leads to definitive quasiparticle 
spectra. We suggest that an elaborate analysis on impurity- induced spectra with quantum interfer- 
ence effect included may be able to pin down the sign and strength of the scattering potential of a 
Zn impurity in low density limit. Numerical simulations calculated with appropriately determined 
impurity parameters are in satisfactory agreement with the observations from scanning tunneling 
microscopy (STM) experiments even in subtle details. 
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Impurity effect has served as a unique probe to the 
mechanism of unconventional superconductivity. Re- 
cently, remarkable improvements in high resolution 
STM experiments provided unprecedentedly delicate 
images of electronic structure around doped Zn in 
BizSrsCaCuaOg+s (BSCCO>i. The defect states in dif- 
ferential conductance map at low temperatures bear on 
almost identical spectroscopic and spatial structures from 
one impurity to the other, with characteristic on-site zero 
bias resonance peak, four-fold symmetry pattern in lo- 
cal density of states (LDOS) and locally suppressed su- 
perconductivity at immediate surroundings of Zn atoms. 
These observations are in qualitative agreement with pre- 
ceding theoretical works which model Zn atoms as iso- 
lated pointlike unitary centers embedded in a supercon- 
ducting background of d-wave like pairing symmetrjiSi^. 
However, although the single-impurity scenario (SI) is 
approximately successful in interpreting the energy and 
spatial symmetry of resonant states, it fails to reproduce 
the subtle features in the spatial profile of LDOS ex- 
tracted by STM measurements. Attempts to reconcile 
this discrepancy include postulations on the tunneling 
matrix of STM probes^, possible Kondo physics from 
staggered moment^ and etc.; nevertheless, the determi- 
nation of the attributes of Zn impurities, with which we 
can substantiate reliable calculations to compare with 
experiments, remains divergent itself: even assuming Zn 
impurities are purely potential scatterers, one could still 
reach completely opposite conclusions that they are ci- 
ther repulsive (as desired to yield unitary limit with re- 
alistic cuprates band) or attractive (relative to the back- 
ground Cu ions) according to their ionic configuration. 
Most recently, an ab-initio calculation based on density 
functional theory concludes that Zn should be repulsive 
centers to electrons^. 

In the meanwhile, the effect of interference between 
the quasiparicle resonances cannot be neglected. Over- 
lapping between disorder-induced Fridel oscillations al- 
ters local landscape in LDOS and populates the low 
energy excitations, which greatly modify the spectro- 
scopic and transport properties of cuprates. There have 
been numerous theoretical studies on many-impurity 
problerr&SiiSiiLiSii^, with emphasis either on the effect 



of impurity network on Fourier transformed power spec- 
trum which can be used to extract the kinematics of pure 
systems or on the asymptotic behavior of Fermi level 
density of quasiparticle states. Despite the creditable 
achievements of these endeavor in interpreting STM data, 
the lingering unjustification of impurity parameters could 
place their arguments on a shaky ground; moreover, none 
of them have given full attention to the potential of com- 
bining interference effect to discuss the impurity spectra 
and identify the impurity characteristics itself. It is then 
the purpose of this paper to argue that through elab- 
orate investigation on quantum interference effect in a 
fully disordered system, we could nail the identities of 
Zn impurity such as its sign and strength and hopefully 
close the debate on this issue to our best. 

We start by writing down the generic many-impurity 
Hamiltonian as, 
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where cl^{cia) creates (annihilates) an electron with spin 
(T at site n and (/>(n) is the impurity strength on the per- 
turbed site. Randomness of impurity distribution makes 
counting relative coordinates information between all the 
impurities a formidable task. Therefore a systematic and 
reliable manipulation of the interference, yet not losing 
the reality and generality, will be helpful. We tacitly sug- 
gest to fix the first impurity at the center (r — 0) of the 
lattice and keep injecting others while requiring that the 
inter-impurity distance of any pair, i.e., Rij = jn — rj| 
has to be greater than a preset cutoff do. This con- 
straint excludes the possibility of cluster formation and 
further guarantees a uniform distribution of disordered 
sites. The central impurity, which is the focus of our dis- 
cussion later, is believed to be the representative of a ar- 
bitrary one surrounded by other disordered sites that are 
randomly distributed in real samples. Therefore, a sin- 
gle parameter, i.e., the concentration, will be the defining 
parameter instead of the relative orientation and distance 
between impurities. The effect of the quantum interfer- 
ence is then parameterized as a function of the concen- 
tration X (or the number of impurities Nimp) in a fixed 
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lattice size. 

We study this problem using the conventional 
Bogolibuv-de Gennes' (BdG) formalism: 
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where Tiy = —tij + {(f>{ri) — fJ,) Sij is the single particle 
Hamiltonian, ty is the hopping integral, ^ is the chemi- 
cal potential. Ay = Y^^it^jl ~ ''ii''jT^ spin-singlet 
bond order parameter, V is the nearest-neighbor attrac- 
tive potential), m" and vV- are the eigenfunctions of the 
BdG equations, and En is the eigen-energy. The self- 
consistent conditions are applied to solve the BdG equa- 
tions: 
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where f{E) ~ l/{e^^ + 1) is the Fermi distribution. 
We use the following parameters throughout this paper: 
(ty) = t = 150 meV, (fy) = t' = -OM, V = It. This 
prepares us a hole-like Fermi surface with the hole doping 
l-uf = 1- EiA4aCia)/^^Ny = 0.15, i.e., optimally 
doped region and a maximum gap magnitude about 45 
meV. The local density of states takes the following ex- 
pression: 
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where the factor of 2 arises from spin degeneracy. The 
summation in Pi{E) is averaged over x My (in our 
paper 20 x 20) wavevectors in the first Brillouin zone. 

Single impurity We start our discussion with a brief 
review on the single impurity model and plot the result- 
ing spectra in Fig.l with the "blocking effect" (which is 
addressed below) included. We unbiasedly picked up a 
repulsive potential (f) = 3.0eV. Because of the Bogoliubov 
symmetry of quasiparticles, the on-site LDTC is expected 
to display two resonance peaks at Vt^ ~ ±0.015i^[see 
Fig.l(a),(b)]. While the spatial scale and symmetry of 
impurity states observed in STM fit well with the pre- 
diction of SI theory, the relative spectral weight distri- 
bution is completely reverted^iS: experimental images at 
f2 = -1.5 meV— displays an on-site maximum spectral 
intensity, a local minimum on its nearest neighbor sites 
and a second maximum on the next nearest sites, as illus- 
trated schematically in Fig. 1(d). This inconsistency gave 



birth to the conjectures regarding the indirect tunneling 
through the insulating BiO layer: the "filter effect" due 
to the coherent tunneling^ and an alternative "blocking 
effect" due to the incoherent tunnelingi^. In our pa- 
per, we will stick to the argument of the latter, in which 
the "local differential tunneling conductance (LDTC)" 
fetched by STM probes is believed to measure acutally 
the averaged LDOS over the four nearest neighbor sites. 
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where e is ztx or zLy. Fig l.(a),(b),(c) are then plot- 
ted with the "blocking effect" included intentionally. In 
Fig. 1(c), the spatial distribution of spectral intensity 
mimics experimental results under this manipulation. 
However, corresponding on-site spectra [Fig l(a)]deviates 
from the result of Pan et al. remarkably: the former 
shows a sharp peak with excess spectral background at 
negative bias and a second resonance at positive side; the 
latter displays a zero bias conductance peak and a sec- 
ond small peak on negative bias; furthermore. Pan et al. 
showed that the Friedel oscillation along the nodal di- 
rection decays surprisingly faster than that of the antin- 
odal direction. This is rather counterintuitive since the 
subgap resonances in d-wavc superconductors overlap 
with the vanishing continuum and hence are virtually 
bounded, with supposedly extended tails along the nodal 
direction due to the vanishing order parameter. Wc argue 
that with the inclusion of quantum interference effect, all 
this discrepancy will be reconciled provided that the im- 
purity parameter is properly settled. 

Many impurity We investigate the quantum interfer- 
ence effect by embedding 10 impurities in an otherwise 
clean lattice, which yields x = 0.41%, close to laboratory 
terms. The lattice size is of A'^j; x TV^ = 49 x 49. Fig. 2 
enumerates the on-site LDTC for different values of im- 
purity strength systematically. Generally, quantum in- 
terference splits the single impurity resonances into mul- 
tiple peaks, as illustrated in Fig. 2 and extracting dis- 
tinctive information from those peaks directly is not easy. 
However, the zero energy residual LDOS p{uj = 0) is par- 
ticularly useful as it evolves systematically with respect 
to the impurity potential. When </> = — 3.0eV, p{uj = 0) 
is considerably large; as (f> increases but remains negative 
(attractive), p{uj = 0) decreases until it is completely de- 
pleted. This tendency can be understood qualitatively in 
the sense that the sample recovers its clean spectra which 
has vanishing p{uj = 0) when less contaminated. When 
(j) becomes positive and increases [repulsive), p{u) = 0) 
aggregates remarkably and multiple peaks merge into a 
"single resonance peak" when (j) is around the value with 
which unitary limit is obtained in single impurity anal- 
ysis. Fig. 2(d) shows such a strong resonance peak at 
negative bias slightly below zero, whose width and height 
are in fair agreement with the result of Pan et al.. The 
likeness between Fig. 1(a) and Fig. 2(d) will plausibly 
suggests the single-impurity physics as dominating mcch- 
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(C) i2=-0.015 



(d) STM measurements 



FIG. 1: (a) LDTC on the impurity-site. Dashed green: pure 
system; solid black: LDTC with the present of single impu- 
rity; (b) LDTC on the nearest neighbor site. Dashed green: 
pure system; solid black: LDTC with the present of single 
impurity; (c) calculation on spatial distribution of LDTC at 
Qo = — 0.015t; impurity site is placed at the center of the 
squared window, (d) Schematic plot of the STM image of 
one impurity-. 



anism and the fact that Zn impurities in the experimen- 
tal conductance maps appear to be isolated entities with 
nearly undisturbed four-fold symmetry in LDTC seem- 
ingly supports this viewpoint. However, it is worthwhile 
to point out that exact width of this zero-bias resonance 
in Fig. 2(d) and STM experimentsi is ~ 10 meV, an or- 
der of magnitude bigger than what single- impurity theory 
expects but agrees well with the impurity bandwidth in 
unitarity low density limit, i.e., 7 = \/niAoEpi^. Hence 
the figurative resemblance between the on-site LDTC's 
of many-impurity and single-impurity must be the conse- 
quence of quantum interference between dilute impurity 
states and the "homogeneous broadening" effect which 
was proposed by Atkinson et ali^ and is not negligible 
in laboratory terms. The possibility of \4>\ being weak 
(less than 0.4 eV) is excluded, as it will not introduce 
any prominent resonant behavior and observable inter- 
ference effect requires inter-impurity distance as close as 
3a, corresponding to an unreasonably large population of 
disordered sites. This is demonstrated in Fig. 2(e) and 
2(f): when = ±0.3 eV, the on-site LDTC of dirty sam- 
ples are essentially coincident with that of single- impurity 
case. Collection of all the facts above and comparison 
with experiments lead us to believe that Zn atoms in 
BSCCO-2212 be repulsive potentials at least in the ze- 
roth order. 

The exact magnitude of Zinc impurity is still open 
to determination, but favorably it should not differ too 
much from the value which approximates single-impurity 
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(a) (|)=-3.0 eV 



(c) (t)=-L2 eV 



(e) (t)=-0.3 eV 



(b) (^=3.0 eV 



(d) 6f=^ .2 eV 




(f) 6f=0.3 eV 
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FIG. 2: The local differential tunneling conductance (LDTC) 
on the central Zn impurity with 10 impurities embedded in 
the lattice. Left panel: negative potentials (a) -3.0 eV (c) - 
1.2 eV (e) -0.3 eV. Right panel: positive potentials (b) 3.0 eV 
(d) 1.2 eV (f) 0.3 eV. Green curves in (e) and (f) are on-site 
LDTC of single-impurity case. 



unitary limit (which is = 1.2eV here) with the cuprates 
band structure since the unitary scattering due to doped 
Zn is verified by the scattering phase shift of 0.497r in 
experimentai. We then adopted 4> = 1.2eV as the in- 
put of calculations and the results are plotted in Fig. 3. 
The spectrum in channel (a) shows a striking agreement 
with the result of Pan et al. with a sharp resonant peak 
forming at 17 = and a second smaller peak on negative 
bias. Apparently, single-impurity scattering fails to in- 
terpret both the asymmetry of the peaks positions and 
the broadened peak widths. While the staggered mag- 
netic interaction may also address the two-peak struc- 
ture, we would rather believe that the negative peak 
should arise collectively from inter-impurity correlations 
and our numerical inspection confirms that it actually 
persists in a fairly wide region of positive 0's; in Fig. 
3(b), the spatial distribution of the resonance around a 
Zn impurity shows the local minima of the four nearest- 
neighbor sites, and the local maxima of the eight 2nd- 
nearest- and 3rd-nearest-neighbor sites clearly forming a 
"box" around this Zn atom; in Fig. 3(c), the normalized 
LDTC G(r) at resonance frequency is plotted as a func- 
tion of distance away from the impurity along the nodal 
and antinodal directions. The strength of the normalized 
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(a)x=0.41% 




(d) x=0.71% 
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FIG. 3: (a) The LDTC Gi{E) on a Zn impurity at r = 
for X = 0.42%. The dashed hne represents the LDTC far 
away from Zn impurities; (b) The spatial profile of LDTC 
maps at i7 = —2.4 meV;(c) The LDTC (normalized to the 
peak value) versus distance from the Zn impurity; (d)(e)(f) 
The LDTC Gi{E) on the Zn impurity at r = for a: = 
0.75%, 1.0%, 1.83%. 

LDTC in the nodal direction is found to decay faster than 
that along the antinodal direction, reinforcing what is ob- 



served in STM experimentsi. All this above agree with 
STM observations remarkably and are confirmed to be 
robust against the change of (f> within a wide range cen- 
tered at 1.2eV; When impurity concentration increases, 
the low energy excitations are further populated and a 
subgap impurity band is established gradually, which will 
suppress superconductivity eventually and is referred to 
as "impurity band" , in analogy to similar phenomena in 
semiconductors. This is illustrated schematically in Fig. 
3(d)(e)(f), where the zero bias impurity resonance is fur- 
ther broadened and finally into multiple-peak structure 
with excessive spectral weight filling up the gap region 
when impurity concentration is several times bigger than 
the experimental value (0.2% — 0.5%). 

Conclusion: We discuss how quantum interference be- 
tween many impurities can produce qualitatively differ- 
ent spectral features when the impurity variables change. 
With a numerical study on the disorder induced local 
spectra (on impurity site), we would like to close the de- 
bate over the identity of Zn impurity by concluding that 
Zn atoms in Bi2Sr2CaCu208-i-5 are repulsive to electrons 
in nature with a strength close to unitary limit provided 
that other internal degrees of freedom was not considered. 
The discrepancy between STM experiments and the re- 
sults of single-impurity analysis can be reconciled satis- 
factorily by taking the quantum interference and insu- 
lating layer blocking effect into consideration simultane- 
ously. Our numerical calculations with impurity param- 
eter determined in this paper match STM experiments 
up to subtle details. We then would emphasize that the 
results obtained in laboratories should actually be inter- 
preted within the frame of collective quantum interfer- 
ence processes rather than the single impurity physics. 
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